Strengthened PT— symmetry with P^P^ 



Miloslav Znojil 
Ustav jaderne fyziky AV CR, 250 68 Rez, Czech Republic 1 



Abstract 

In Quantum Mechanics working with non-Hermitian PT— symmetric Hamiltonians 
(i.e., with an indefinite metric P in Hilbert space) we propose to relax the usual 
constraint P = P'. We show that this merely induces certain "hidden" symmetries 
responsible, say, for the degeneracy of levels. Using a triplet of the coupled square 
wells for illustration we show that the bound states may remain stable in a large 
domain of couplings. 
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1 VT— symmetric Quantum Mechanics and its two 
alternative scenarios 



VT— symmetric Quantum Mechanics of C. Bender et al [1] studies non-Hermitian 
Hamiltonians H ^ W with the peculiar property 

VT H = HVT . (1) 

In the original and simplest one-dimensional version of the theory [2] the symbol 
T denotes the complex conjugation, i.e., an antilinear involution with the prop- 
erty TidtT = —idt interpreted as time reversal. The Hamiltonians themselves 
are assumed symmetric so that we may replace H — > h = h T and deduce that 
ThT = h* = h), T hVT = WV* and 

h) = V^hV = V*h (V*)" 1 . (2) 

The choice of the parity V with the properties V = V* = V^ 1 = V T = V^ in 
ref. [2] inspired A. Mostafazadeh [3] who recommended a transition from the physics- 
inspired symmetry (1) to its mathematically equivalent representation in the form 
of the V— pseudo-Hermiticity requirement 

H ] =VHV~ l . (3) 

He emphasized that eq. (3) may be read as an isospectrality property where one might 
work with asymmetric Hamiltonians and with the "generalized parity" operators 
which need not be involutive at all, P-^P^ P _1 . 

In such a perspective the boldface symbol P may represent an arbitrary auxiliary 
operator. The involutive character of the Hermitian conjugation in eq. (3) implies 
that 

H = (tft) f = (pt)- 1 #tpt = (pt)- 1 Ptfp-ipt ( 4 ) 
which gives us the two alternative possibilities; 

[a] a self-adjoint pseudo-metric P = is chosen, or 

[b] non-Hermitian operators P ^ P' are admitted. 

In the light of the current literature the "trivial symmetry" scenario [a] seems to be 
"the only useful" option where P becomes a pseudo-metric (often called "indefinite 
metric") in the physical Hilbert space of states H. 
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We believe that the non- Hermit ian alternative [b] may prove equally interesting. 
Indeed, there is no real reason for ignoring the class of operators S = P" 1 7^ / 
which represent a nontrivial symmetry (4) of the Hamiltonian. In what follows we 
intend to support such a point of view by an illustrative construction. For this 
purpose we shall pick up one of the most elementary coupled-channel symmetries 
S ^ I derived from the non-involutive and non-Hermitian toy operator 



Provided that its sub-operator V remains defined as the parity, Vip(x) = ip(—x), it 
cannot be interpreted as a metric because its eigenvalues are complex. 

A supplementary reason for the present use of the non-Hermitian P of eq. (5) 
emerges once we return back to the symmetric Hamiltonians h = h T in eq. (2) 
where any "early generalization" V — > TZ of the parity would lead immediately to an 
alternative explicit constraint 



We see that the "alternative generalized parities" 1Z would have to be unitary when- 
ever one assumes that they are not asymmetric. One should keep in mind that the 
latter consistency constraint definitely differs from its predecessor eq. (4). Of course, 
it can again be interpreted as imposing an additional symmetry upon the Hamilto- 
nian. Thus, one feels that eq. (3) with self-adjoint V need not necessarily offer the 
only productive way towards a generalization. 

2 Non-Hermitian triplet of coupled square wells 

Equation (5) and the non-Hermitian Hamiltonian of the triple-channel form 




/ V \ 




\ V J 



TZTZ*h = hWl*. 



(6) 



H = 



U a (x) 

V V a (x) 



dx 2 



U b (x) 
3 + D b (x) 
W a (x) 



dx 2 



V b (x) \ 
W b (x) 

3 + D c (x) ) 



(7) 



will be assumed inter-related by our present modification 



H ] = PHP \ 



P ^ P f 



(8) 
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of the P— pseudo-Hermit icity condition (3) re-written in the form adapted to the less 
common scenario [b]. Our choice of the interaction potentials will be dictated by 
the exact solvability requirement in a way inspired by the simplicity of the various 
single- channel square- well models [4] - [6] in scenario [a]. 

No innovations will occur in the real, Hermitian part of the present potentials, 

ReD aAc (x) = ReU a , b (x) = ReV a , b (x) = ReW a , b (x) =0, x G (-1, 1), 
ReD aAc (x) = ReU a , b (x) = ReV a ,b(x) = ReW a , b {x) = oo, x (-1, 1). 

In the same spirit, the imaginary potentials in Hamiltonian (7) will be postulated 
piecewise constant. Their specification 

lmU a , b (x) = lmV a!b (x) = lmW a>b {x) = Y>0, x G (-1,0), 

lmU a , b (x) = lmV a , b (x) = lmW a , b (x) — — Y, x G (0, 1), (10) 

lmD aM (x) = Z, xG(-l,0), Im D aAc (x) = -Z, x G (0, 1) 

containing two free real coupling constants results from the pseudo-Hermiticity con- 
dition (8) and from our choice of the generalized parity (5). This defines the family 
of the coupled-channel Schrodinger equations 
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) 



where the energies E may be assumed real, for the small couplings Y and Z at least 
[7]. As long as we put our system in a box (9), the elementary boundary condition 









( ) 


if b (x) 









V <Pc(x) J 




x=±l 


VO ) 



determines all the bound states of the model. 



3 Solutions and the degeneracy of their spectrum 

The obvious ansatz for the general solution 



. Ci a ' M sin KL (x + l), XG(-1,0), 

VW^) = <^ _(„ )6)C) . , c fn ^ (13) 

C R ' sin kr{— x + 1), x G (0, 1) 
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takes into account the boundary conditions and necessitates only the appropriate 
matching of all the three wave functions (p a ,b,c( x ) and of their first derivatives at 
x — 0. This imposes the six complex constraints 



/~i(a,b,c) ■ /-y(a,b,c) 

C L sin kl = Cr sin k r , 



K L Cl' ' COS K L = -Kr Cr ' COS Kr. 



The first triplet merely defines, say, constants C^"' 6 '^ as functions of Cr'^ and kl,r- 
The ratio of both of these equations eliminates all these constants and leads to the 
single complex equation 



k l cotan k l = — kr cotan kr. 



(14) 



Under this constraint our final quantization condition will result from the insertion 
of the ansatz (13) in the Schrodinger eq. (11) at x E (—1,0), 



(15) 
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and at x E (0, 1), 
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(16) 



As long as the energies are assumed real, this indicates that we may put kr — s+it — 
k* l with, say, positive s > and any real t E (— oo, oo). In this notation the complex 
eq. (14) degenerates to the real implicit formula 



2s sin 2s + 2t sinh 2t = 



(17) 



which first occurred in ref. [8] and which has thoroughly been studied in ref. [9]. In 
our present model, eq. (17) has to be combined with the complex secular equation 



det 



( (s + itf-iZ-E -iY 

-iY (s + it) 2 -iZ-E 

\ -iY -iY 



-iY \ 
-iY 

(s + it) 2 -iZ-E ) 



0. (18) 



The resulting three real conditions have to specify the two unknown real parameters 
s — s n , t — t n and the energy E = E n , n = 0, 1, . . .. Once we set E 



s 2 -t 2 - a 



4 



(with a real a) and Z = 2st + (5 (with a real (5) we may re-write eq. (18) as a pair 
of the real polynomial equations 

a 3 -3a ((3 2 -Y 2 ) = 0, , s 

VP ; ' 19 

/? 3 -3/3 (a 2 + F 2 ) + 2F 2 = 0. 

In the preliminary test we shall assume that a = (X(tentative) 0- From the first row 
we get af tentative ) = 3 (f3 2 — Y 2 ) which simplifies the second row to the solvable cubic 
equation with the three roots, 

^(tentative) _ v ^(tentative) _ n(tentative) _ 1 v (tentative) / n 

Pi — 1 1 l J 2 — H3 ~ ~2 

Their insertion in the definition of a gives the respective solutions 

3i 

ai = 0, a 2 ,3 = ±—Y 

all of which contradict our assumptions. We are forced to conclude that we always 
have the vanishing a — 0, 

E = s 2 - t 2 . (20) 
At a = the secular equation (19) leads to the unique triplet of roots 

A = -2Y, (3 2 = fo = Y, a = 0. (21) 

Their respective insertion in eqs. (15) and/or (16) gives the unnormalized eigenvector 

(cg,C$,cg)~ (1,1,1) (22) 

plus the two other, due to the degeneracy, non-unique vectors available, say, in an 
orthogonalized representation 

(Cg, Cl% Cg) ~ (1, -1,0), (^J, Cg eg) ~ (1, 1, -2) (23) 

which, incidentally, coincides with the Jacobi-coordinate recipe for the three equal- 
mass particles [10]. 

Our knowledge of the roots (21) enables us to eliminate one of the real parameters 
(say, t) as lying on one of the two different hyperbolic curves, 

t = &\s) = ±-Z eff (a), Z eff (l) = Z + 2Y, Z eff (2,3) = Z-Y. (24) 
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Our construction of bound states is completed. In terms of the parameters s and t 
they are determined by formulae (13), (20), (22) and (23). The parameters them- 
selves must be fixed by the pair of eqs. (24) and (17). In a way described more thor- 
oughly in ref. [8], all the physical roots of our secular eq. (17) may be re-parametrized 
by the formula 

8 = S n = ^llpl + (_l)n eB> n = ,l,... (25) 

where the new unknown quantities e n remain small not only in the vicinity of the well 
known Hermitian case where both the coupling constants Y and Z remain sufficiently 
small but also at all the sufficiently large n > no- Thus, one may calculate them 
perturbatively in both these regimes [11]. 

4 The determination of the domain where the en- 
ergies remain real 

In the s — t plane we may visualize all the roots (s n ,t n ) as intersections of the 
two hyperbolic curves (24) with all the family of the (t — > —t) — symmetric ovals 
representing the complete graphical solution of our second secular implicit constraint 
(17) (cf. [8]). As long as Y > 0, all the present energies E n remain real if and only if 

Y - Z crit <Z< Z crit - 2Y. (26) 

We may recollect the commentary in ref. [9] and conclude that Z crit 4.475 at 
the present choice of the units h = 2m = 1. In the (Y, Z) plane the condition 
(26) is satisfied inside a fairly large triangle with (approximate) vertices (0,4.4753), 
(0,-4.4753) and (2.9835,-1.4918) (cf. Figure 1). 

The critical parameter Z crit remains the same for several different square-well 
systems. It determines the boundary of the physical domain in the single-channel 
square well as well as in all its classical [11] and supersymmetric [12] partners. Still, 
only its four-digit estimate has been published up to now [9]. Moreover, even that 
improvement of the original three-digit estimate of ref. [8] by one digit did not 
seem easy. This apparently discouraged, undeservedly, any other attempts. For 
example, an interesting alternative approach of ref. [13] using a discretization of the 
coordinates has only been studied in the lowest possible approximation giving just a 
schematic initial estimate 4^/2 ps 5.66 of Z crit . In fact, a more complicated problem 
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with complex energies has been solved during the most successful numerical attempt 
in ref. [9]. For this reason, let us show now that a systematic improvement of the 
precision of Z crit can be made feasible at a reasonable computational cost. 

Firstly let us summarize the situation where one locates, graphically [8] , the first 
two single-channel non-Hermitian square-well energies £0,1 = s§ x — ^0,1 as related 
to the neighboring intersections (so,ij*o,i) °f the implicitly defined Z— independent 
curve s = s (a) (t) [with 2s (a) (t) sin[2s (a) (t)] = -2t sinh2t from eq. (17) above] with 
the Z— dependent but much more elementary hyperbolic curve s^(t) = Z/(2t) of 
eq. (24). In such a graphical setting it was clarified that the maximal Z = Z crit at 
which both the energies £0,1 remain real is the point at which they [as well as the 
neighboring intersections (so,i, £0,1)] merge. Thus, the value of Z = Z cr n is defined 
as a parameter of confluence at which Eq precisely coincides with E\. 

At this point the curves s^(t) = n — e(t) and s^(t) will osculate at a certain 
"intersection" point t crit and "exceptional" energy E crit . We have to guarantee that 
both our curves and both their tangents coincide, 

Z ■ Z 
e(tcrit) = 7T - — — , d t e(t crit ) = — f^. (27) 

A* crit Al cr n 

The derivative is defined in terms of the positive shift function e(t) < 7r/4, 

sinh 2t + 2t cosh 2t 
tS ^' ~ 2 [tt -e(t)] cos2 £ (t) -sin2e(i) 

but the definition of the function e(t) itself is merely implicit, 

t sinh 2t 

sm [2e(t )] = — -. 

7r — e(t) 

Fortunately, it may be re-interpreted as a quickly convergent iterative recipe of a 
'generalized continued-fraction' type, 



, N 1 
Zinewft) = arcsm 



t sinh 2t 



[28) 



With the initial £(i ower )(t) = 7r/4 and £( Mpper )(t) = 0, the N—th iteration of eq. (28) 
represents the desired explicit definition of the respective approximate functions 
£(i ower ){t) and £( upper )(t). Their knowledge enables us to solve the two coupled alge- 
braic equations (27) for the two unknown quantities t cr u and Z crit at each particular 
choice of N. The both-sided convergence of this recipe is illustrated in Table 1. For 
the sake of completeness, its last-row items may be complemented by the correspond- 
ing values of t ( ^ er) = 0.839393459, t^ per) = 0.839393461, s^ er) = 2.665799044, 
s ^ er ) = 2.665799069 and E«™ er) = 6.401903165 and E^ er) = 6.401903294. 
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5 A remark on the interpretation of the model 



In summary, we felt inspired by several physical applications of scenario [a] which 
have recently been offered within relativistic quantum field theory [14] and first- 
quantized relativistic quantum mechanics [15] as well as in quantum cosmology [16] 
and in the classical magnetohydrodynamics [17]. In these cases one often employs the 
partitioned and manifestly Hermitian and involutive P, in the latter three contexts 
at least [18]. In our present letter we complemented these studies by an illustration 
of several new possibilities emerging within the scenario [b]. 

In our present coupled-channel model the Hilbert space is partitioned in sub- 
spaces. In fact, there is no real novelty in such a procedure of the model-building 
as the various P— pseudo-Hermititian coupled-channel operators are known to re- 
sult from the relativistic Sakata-Taketani equations [19, 20] and from their various 
higher-spin generalizations and/or non-relativistic analogues [21]. A fresh example 
of a coupling of channels in non-Hermitian context may be found in our recent Klein- 
Gordon study of the influence of certain external solvable delta-function interactions 
[22]. The Hermitian part of our forces was postulated there in the same deep square- 
well form of eq. (9) and only the pseudo-metric P = has been chosen Hermitian, 
of type [a] . 

Let us emphasize that from the point of view of Quantum Mechanics it is im- 
portant to know that there exists (at least one) specification of the scalar product 
which leads to the positive definite physical norm. Unfortunately, its explicit con- 
struction is usually fairly complicated in practice. In this sense, the exact solvability 
of the square-well-like models simplifies significantly the perturbative construction of 
the related "physical" metric 6 ^ P in Hilbert space [11]. In the other words, our 
knowledge of G enables us to assert that all the observables in our models acquire the 
so called quasi- Hermiticity property [23] and that only in terms of the (by definition, 
positive-definite and self-adjoint) they may be assigned the standard probabilistic 
interpretation [23] - [28]. 

In our present non-relativistic coupled-channel example the transition to the sce- 
nario [b] did not influence the methods of the construction of the physical metric so 
that they need not be discussed separately. Even all of their technical details remain 
the same for our specific choice of the operator P since P 3 = V. This is an accidental 
aspect of our assumptions (5) and (8) which makes the validity of the standard rule 
(3) an immediate consequence of our assumption (8). 
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In an opposite direction one may say that our present assumptions concerning the 
symmetry of the Hamiltonian are stronger since our H commutes with the product 
S = [P _1 ]*P 7^ /. In effect we preserve and complement the "old" eq. (3) by another 
assumption. As long as we impose more symmetry, the degeneracy of some levels 
occurs. From the practical point of view such a phenomenon is the consequence of 
our choice of the complicated P ^ Pt Of course, a wealth of new features of the 
spectrum may be expected to emerge from the more systematic study of some less 
schematic non-self-adjoint "generalized parities" P in the future. 
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Table 1: Numerical determination of the critical coupling 



iteration 
N 


ry (lower} 


r? ( uwer) 
^crit 





4.299 


4.663 


2 


4.4614 


4.4857 


4 


4.47431 


4.47601 


6 


4.475239 


4.475357 


8 


4.47530381 


4.4753119 


10 


4.475308262 


4.475308823 


12 


4.475308560 


4.475308614 



Table captions 

Table 1. Numerical determination of the critical coupling 
Figure captions 

Figure 1. Triangular domain of the allowed couplings Y and 
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